Thermoelectric properties of junctions between metal and models of strongly 

correlated semiconductors 



Massimo Rontani 1 ' 2 '^ and L. J. Sham 1 

1 Department of Physics, University of California San Diego, La Jolla, California 92093-0319 
2 INFM National Research Center on nano Structures and bioSystems at Surfaces (S3) and 
Dipartimento di Fisica, Universita degli Studi di Modena e Reggio Emilia, J^llOO Modena, Italy 

(Dated: February 2, 2008) 
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potential of the junction as a thermoelectric device is discussed. 
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I. INTRODUCTION 

New thermoelectric coolers and power generators are 
under massive investigation, i These devices are in general 
even more reliable than commercial heat-exchange refrig- 
erators but their efficiency, in the best cases, is much 
lower. The quality of the material need in a thermo- 
electric device is defined by the dimensionless figure of 
merit ZT, where T is the absolute temperature, and Z 
is expressed in Eq. Q in terms of transport coefficients. 
Currently, the highest value of ZT, which is ~ f at room 
temperature, is found in Bi-Te alloys^ New devices would 
be competitive with traditional refrigerators if ZT were 
about 3^4: the major lack of high-ZT materials is at 
temperatures below 300 K. 

In this paper— we propose a junction of metal and 
strongly correlated semiconductor as the basis for a possi- 
ble efficient low-temperature thermoelectric device. This 
system embodies previous intuitions that were recog- 
nized fruitful 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 : 14 ! 15 ! 16 in different materi- 
als/devices such as rare-earth compounds, superlattices, 
and metal/superconductor junctions. We now outline 
these concepts and briefly review the related literature. 

From the definition 

(K e + Hi) 

where Q is the absolute thermopower (Seebeck coeffi- 
cient), a the electrical conductivity, and n e and ki the 
electronic and lattice part, respectively, of the thermal 
conductivity, it follows that an ideal thermoelectric ma- 
terial should have high thermopower, high electrical con- 
ductivity, and low thermal conductivity. Semiconductors 
seemed to have the optimum collection of these proper- 
ties, in contrast with metals, which have high a but low 
Q, and insulators, which have high Q but low a. How- 
ever, in the last thirty years no substantial enhancement 
of ZT beyond ~ 1 was obtained. 

A breakthrough came with the synthesis of new mate- 
rials such as filled skutterudites, which have nearly bound 



rare-earth atoms, closed in an atomic cage, whose "rat- 
tling" under thermal excitation scatters phonons, then 
dramatically reducing More generally, heavy atoms 
in compounds help with lowering ki. Mahan and Sofc 5 
also showed that the best bulk band structure for high Q 
is one with a sharp singularity in the density of states very 
close to the Fermi energy. These results provide the first 
idea in the search for the best thermoelectric, namely 
to look at rare-earth compounds as major candidates. 
In fact, mixed valence metallic compounds (e.g. CePd3, 
YbAl 3 ) show high values of Q, but at the present time 
no useful value of ZT has been reported^ 

The second idea is that the best thermoelectric must 
have an energy gap. Because in ordinary semiconductors 
the optimum band gap is predicted to be about 10 fee T 
(fee being the Boltzmann constant), 6 one is led to con- 
sider small-gap semiconductors for low-temperature ap- 
plications. If the chemical composition of semiconductors 
includes transition metals or rare earths, conduction and 
valence bands are frequently strongly renormalized by 
correlation effects, forming a temperature-dependent gap 
(see Fig. |2Jl : this is the case for mixed- valent semiconduc- 
tors, usually cubic, whose relevant electronic properties 
may be modeled by a /-flat band and a broad conduc- 
tion band, with two electrons per unit cell. 18 This class 
of materials consists of two subclasses: The first is the 
Kondo insulator, characterized by very strong Coulomb 
interaction between electrons on the same rare-earth site, 
usually described by the slave boson solution of the An- 
derson lattice hamiltonianii 9 * 2 ^ Mao and Bedell predicted 
a high value of ZT for bulk Kondo insulators (the lower 
the dimension, the higher the value) £ however, some ex- 
perimental reports seem to exclude this possibilityiSii 2 ^ 
The second, called the electronic ferroelectric (FE)j2i 
consists of semiconductors with high dielectric constants, 
such as SmBg and Sni2Se3, and it is modeled by the 
self-consistent mean-field (MF) solution of the Falicov- 
Kimball Hamiltoniani^i The ground state of the insu- 
lating phase is found to be a coherent condensate of d- 
electron//-hole pairs, giving a net built-in macroscopic 
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polarization which breaks the crystal inversion symmetry 
and makes the material ferroelectric ^ 

Another useful observation is that ZT is reasoned to 
increase in quantum-well superlattices, due to the mod- 
ification of the density of states&tLiSi Moreover, super- 
lattices with large thermal impedance mismatch between 
layers seem very efficient at reducing ni because interfaces 
scatter phonons very effectively^ However, if the trans- 
port is parallel to the layers, the parasitic m of the bigger 
gap materialise or the tunnelling between conduction 
layers 2 ^ can drastically decrease ZT. Besides, if intra- 
layer transport is diffusive and described by bulk param- 
eters, ZT for the whole superlattice cannot be higher 
than the maximum value for the single constituents^ 

In addition to the above literature, this work was 
stimulated by some recent advances in thermoelectric 
applications of junctions. Nahum and coworkers^ 
built an electronic microrefrigerator based on a metal- 
insulator-superconductor (NIS) junction. Subsequent 
experimental 28 and theoretical workiiiSLSI confirmed 
this new idea. Edwards and coworkers^ showed that 
tunneling through structures with sharp energy features 
in the density of states, like quantum dots and NIS junc- 
tions, can be used for cryogenic cooling. It seems then a 
natural extension to us to study the junction between a 
metal and a strongly correlated semiconductor. 

There were recent proposals for devices such as semi- 
conductor/metal superlattices with transport perpendic- 
ular to interfaces. Mahan and WoodsM suggested a mul- 
tilayer geometry with the thickness of the metallic layer 
smaller than the electronic mean free path and the semi- 
conductor acting as a potential barrier (thermionic re- 
frigeration). Independently, Moyzhes and Nemchinsky 15 
proposed a similar configuration with the metallic layer 
thickness comparable with the energy relaxation length. 
We cite also Min and Rowe's ideai& of using Fermi-gas / 
liquid interfaces. While the analysis of electronic trans- 
port of Ref. ^3 an d \IM is n °t applicable to our study, 
because of the sharp energy profile of the transmission 
coefficient across the junction which we examine, and of 
the nature of the correlated semiconductor ground state, 
this experimental geometry, with the strongly correlated 
semiconductor replacing the barrier layer, could be im- 
plemented, as we suggest at the end of Sec. IIVCI 

In sum, we build on two key ideas from the above lit- 
erature for increasing ZT. One is to utilize the sharp 
energy features in the density of states of bulk materi- 
als as in strongly correlated semiconductors. The other 
is to exploit the good thermoelectric characteristics of a 
junction. In this paper, we combine these ideas in ex- 
ploring the thermopower behavior of a junction between 
a metal and different classes of semiconductors. We de- 
scribe the gapped material on one side of the junction 
as the solution of the Falicov-Kimball HamiltonianSi in 
different regimes. In particular, we consider the case of: 
(i) Electronic Ferroelectric (FE), where, because of the 
Coulomb interaction between /-holes and d-electrons, the 
MF insulating ground state is a condensate of excitons; 23 



(ii) Narrow Band semiconductor (NB), characterized by 
the d-f band hybridization: it would be a Kondo insu- 
lator had we take into account the /-/ electron repul- 
sion; (iii) Broad band semiconductor (SC), for compar- 
ison. In these three cases, we solve the electronic mo- 
tion across the junction by means of a two-band model 
analogous to the Bogoliubov-de Gennes equations^ in a 
finite-difference form. In addition to the clean interface, 
we consider an "impurity" overlayer made either of rare- 
earth atoms, with relevant atomic orbitals of /-type, or 
of atoms with d-type orbitals, like transition metals. In 
the latter case, electrons can hop from these "d-impurity" 
sites to adjacent neighbor atoms, while in the former /- 
impurity case hopping is assumed negligible. This sce- 
nario is motivated by recent advances in atomic layer 
fabrication. We compute ZT for the interface via a lin- 
ear response. We find that ZT can be greatly enhanced 
by the presence of a suitable /-impurity layer at the in- 
terface. In particular, for a fixed working temperature of 
the junction, an optimum energy of the /-impurity level 
exists which maximizes ZT, especially at low tempera- 
tures. In these regimes, bulk thermal conductivity would 
be dominated by phonons which would reduce ZT. How- 
ever, one can fabricate the junction with two materials 
with large thermal impedance mismatch, so that phonon 
scatterings at the interface decrease the thermal conduc- 
tivity. Thus, phonon conductivity would not diminish the 
high ZT found. To this aim we propose an experimen- 
tal setup, namely perpendicular transport in metal/FE 
superlattice. 

The structure of the paper is as follows: in Section 
UTI we describe the model Hamiltonian, in Sec. IIII Al we 
solve the electronic motion across the junction, and in 
Sec. IIII 51 we compute transport coefficients and ZT. In 
Sec. IIVI we present and discuss our results, for the clean 
interface (|I V A\i . and for the d- (|IVB|) and /- (jIVC|) im- 
purity layer. Also we briefly discuss some structural re- 
laxation effects of the interface fsec. HVDll . Our conclu- 
sions are in Sec. Ivl 

II. THE MODEL 

We introduce the one-dimensional spinless Hamilto- 
nian Ti. to model the motion across the junction along 
the z direction perpendicular to the interface between a 
metal and different types of semiconductor. TL is given 
by the sum of three terms: 

W = W mc tal + ^interface + ^FK- (2) 

Wmotai is a tight-binding Hamiltonian describing the 
metal on the left side of the junction: 

H mcta i = (e' d + eV) £ d)dj -t'J2 d ] d ^ + H - c - ( 3 ) 

j<0 j<0 

Here dj destroys an electron of charge e < at the lat- 
tice site with energy s' d and position z — aj, a is the 
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FIG. 1: Pictorial representation of the junction along the z- 
axis perpendicular to the interface: the atomic sites at z — aj 
have d-type orbitals while those at z = a(j + 1/2) have /-type 
orbitals, even and odd under spatial inversion, respectively. 
Different symbols for the impurities at z = and z = a/2 are 
drawn; also the parameters of the Hamiltonian of Eq. are 
depicted. 



lattice constant, t' is the hopping parameter for nearest- 
neighbor sites, and V is the electrostatic potential across 
the junction by applying an external bias. In our model, 
V is constant in the two bulk regions with a discontin- 
uous step at the interface, although the actual profile 
should be determined self-consistently together with the 
electronic density. 

The third term Hfk m Eq. J5J) is the Falicov-Kimball 
Hamiltonian24 referring to the insulator on the right side 
of the interface. In addition to the Coulomb interac- 
tion Udf between the d-electrons at sites z = aj and 
the /-electrons localized on atoms at z = a(j + 1/2), 
we add a hybridization term (Vdf is the modulus of the 
hybridization integral) between d and / orbitals. Be- 
cause the crystal has inversion symmetry, this term has 
to be odd, namely the hybridization term between sites 
at z = ja and z = (J + 1/2) a and that between sites at 
z = (j + 1/2) a and z = (j + 1) a have opposit sign (see 
Fig.[I): 



Hfk = Hfk 



n 



(2) 
FK' 



(4) 



n 



(i) 

FK 



E + £ ~/ E /J+l/a/3+1/2 - 1 E 4 d 3 + 1 

j>0 j>0 j>0 

v v E f]+i/2 d i + v <» E + H - c - ( 5 ) 

j>0 j>l 



KfK - U df E djdjfj +1 /zfj+l/2 
j>0 

+ ^/E4^/]-i /2 /i-i/2, (6) 

J>1 



where t is the hopping coefficient, and if and id are the 
/- and d-site energies, respectively. 

The second term ^interface in Eq. J2J is the Hamilto- 
nian at the interface, which describes the overlayer made 
of d-sites at z = and /-sites at z = a/2: 



H 



i:icrfacc — ^dodftdo + ifi/2fy 2 fl/2 

+ U d f L dld fl /2 h/ 



UdfRd\diflj 2 f 1 /2 
-V d f L fl /2 d + Vdfufl^d! + H.c. (7) 

Here we have included the possibility of "impurity" 
atoms at the interface, namely one with average energy 
ido at position z = and another one with energy inn 
at position z = a/2. Since the impurity atoms have or- 
bitals differing from those in the bulk, the associated Udf 
and Vdf parameters generally change. We denote them 
by U d fL, VdfL, U d fR, and V d fR, referring to the couples 
of sites z = 0, a/2 and z = a/2, a, respectively (Fig-QJ. 

We give the MF solution of the Hamiltonian H of 
Eq. (J2J, by means of an approach (see Sec. 1111 Aj l anal- 
ogous to the Bogoliubov-de Genncs method 31 (BdG). In 
contrast to the electron-electron pairing in the Bardeen- 
Cooper-Schrieffer (BCS) MF theory of superconductiv- 
ity, the pairing here occurs between d-electrons and /- 
holes. To proceed, we assume that 



dldj) =n d Vj > 1, 



Udf (d)f j+1/2) = A Vj > 1, 
Udf (4/,--i/a) = A Vj>2, 



4 rf o) = n d0 , (/1/2A/1 



1/1/2, 



= A/ 



(8) 



(9) 



(10) 



fill 



(12) 



(13) 



Here (. . .) is the symbol for the quantum statistical av- 
erage. Note that, in addition to the usual mean or- 
bital occupations rid and n/ of standard Hartree-Fock 
theory (0 < n d ,nf < 1), we also introduce the non- 
vanishing pairing potential A, characteristic built-in co- 
herence of the d-electron//-hole condensate. Follow- 
ing de GenneSfi^i we now write an effective Hamiltonian 
H MF , to be computed self-consistently together with the 
energy spectrum: 



H 



MF 



H-n 



ctal 



ajMF 

n 



interface 



-i/MF 
n FK > 



(14) 
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ft 



MF 

interface 



e do d f d + e f 1/2 fl /2 f 1/2 - &Lfl/ 2 do 
Ai?/| /2 di - V dfL fl /2 d Q + V d f R f\ /2 di 



-UdfLndO n fl/2 



A, 



H.c. 

(15) 



n 



MF 
bulk 



= Sdid\di + e d d]dj + e f ^ /J+i/ 2 /j+i/2 



j>Q 



i>0 



d]d j+1 - V df J2f j+1/2 dj 



J>0 



j>l j>Q 

A^/ j f _ 1/2 d J + H.c. - U dfR n d n fl / 2 



U df n d n f 



j>1 



2U df n d n f 



E 2 i A i 2 ' 

i>o 



(16) 



Here we have defined the renormalized energies 

£<iO = £dO + U d fLTLfi/2, 

e /i/2 = £/i/2 + U d fLn d0 + U d f R n d , 

£dl = id + U d fR fljx/2 + U d f Uf, 

s f = if + 2U df n d , 
s d = i d + 2U d fn f . 



(17) 



We shall treat the renormalized quantities e d o, £/i/2j £di, 
Sf, and e d as material parameters. For simplicity, we 
shall assume t' > t, e d \ = e d , and e' d — e d — 0, that is 
the middle of the d-band on both sides of the junction is 
aligned, but the metal bandwidth (4t') is larger than the 
semiconductor bandwidth (4i). We will consider only 
the case e/ = 0, that is the fiat band is in the middle 
of the d-band in the right-side material (see Fig. |3J). 32 
With these assumptions, and dropping all the constant 
terms in Eqs. (|15I16|) . we obtain the expression of the 
Hamiltonian 7i MF which we will use in the following: 



n 



MF 



5V}^d]dj + £ d od f Q d Q + e fl / 2 fl/ 2 f 1/2 

j<0 

- t'J2d]d j+ i-tJ2d]d j+ i 

j<0 j>0 

-(Az, + V dfL )fl /2 d Q - (A* - V dfR )fl /2 di 
-(A + ^)^/j +1/2 ^ 

j>0 

-( A -%)Eii^+ ac - 

Note the different parity of A- and Vd/-terms: while the 
Vd/-term is odd, the A-term is even^S 

The values of averages on the left side of equations 
(|8I13|) are spatially inhomogeneous and should be deter- 
mined self-consistently and simultaneously with the solu- 
tion of the BdG equations. We make the approximation 
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FIG. 2: Quasi-particle bulk energy bands of materials on 
both sides of the junction vs wave vector at T = 0, for the 
three cases under study, in the "semiconductor representa- 
tion." All energies are in units of t, and e = corresponds to 
the chemical potential \i. The second panel [case (i)] shows 
also the "bare" band structure (dashed lines), i.e. that not 
renormalized by correlation, with Vdf = A = 0. Since the 
alignment of fj, with respect to the left-hand side metal band 
depends on the kind of bulk semiconductor on the right-hand 
side, here we assume a metal/FE setup, with metal (semicon- 
ductor) bare bandwidth 4t' (4i), and t' = 5t. Note that in the 
case (i) [electronic ferroelectric] and (ii) [narrow band semi- 
conductor] the gap is indirect, while it is direct in the case (iii) 
[broad band semiconductor] . Besides, in case (i) [(ii)] the bot- 
tom (top) of conduction (valence) band is much flatter than 
the top (bottom) of valence (conduction) band, while in case 
(iii) the curvature of the two bands close to the direct gap is 
comparable. The trasformation k — > iv/a — k, e —* — e maps 
the FE electron (hole) band into the NB hole (electron) band. 
In this plot parameters are: (i) A(T = 0)=0.49t, V df =0. (ii) 
A=0, V df =0A9t. (iii) t"=0.141t, K=0.071t. 



that those are the bulk averages appropriate to each side 
of the interface layer. In particular, we take the order 
parameter A as constant on the right side of the inter- 
face, and zero on the left side^ The quantities A^, A R 
are used to characterize the effect of the interface on A. 

The temperature dependence of A is much stronger 
than that of the average occupation of orbitals. We de- 
termine numerically A, together with the chemical po- 
tential /i, as the bulk value deep inside the right material 
in a self-consistent way. For this purpose we use the bulk 
Hamiltonian H^ k : 



ivMF 
rl bulk 



-t J2 d)d j+1 - (A + V df ) f} +1/2 d, 



(A-V d f)^2fl 1/2 d j+ H, 



(19) 



Here the index j runs over the whole space. We always 
assume an electronic occupation of one electron per unit 
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cell (a cell contains one ci-site and one /-site), namely 



A. The Bogoliubov-de Gennes equations 



1 



3 



d)d 3 - 



+ f j+ y 2 fj+l/2 



= 1, 



(20) 



where N s is the number of cells. Eq. I|2U|) allows us to cal- 
culate the chemical potential /i for the bulk. Because the 
effective Hamiltonian of Eq. (|19(l is strongly temperature- 
dependent, so /x is. Since /i refers to the bulk value inside 
the semiconductor, the metal on the left side acts as an 
electron reservoir. 

We will consider three specific cases for the actual val- 
ues of parameters in the Hamiltonian H MF of Eq. (|18J) : 
(i) V d f -» and non-zero U df (V dfL = V dfR = 0), 
i.e. the case FE of the electronic ferroelectric with neg- 
ligible hybridization. 34 (ii) Udf — + and non-zero V d f 
(U d fL — U d fR = 0), i.e. the case NB of a /-band hy- 
bridized with a (i-band, that is a narrow band semicon- 
ductor which would be a Kondo insulator if it were driven 
by the condensation of slave bosons, (iii) A ordinary, 
"broad band" semiconductor SC, without center of in- 
version. In the latter case we still use Tt MF but we re- 
gard it simply as a one particle Hamiltonian where the 
"d" and "/" indices are pure labels, and we rename V d f 
as t, considering it as the odd part of a hopping coeffi- 
cient between nearest neighbor sites, A as 14, the even 
part, and t as t" , a second nearest neighbor hopping co- 
efficient. We make the choice t ^> V e , V e ~ t", leading to 
broad conduction and valence bands with different effec- 
tive masses. 

While in case (i) A has to be determined self- 
consistently and in case (ii) Vdf is a material parameter, 
there is a formal canonical transformation connecting (i) 
to (ii), 



d 3 ^d]c mj , f j+1/2 -> /T +1/2 c 

eV -> -eV, 



-iTrO' + l) 



E/l/2 -* ~ £ fl/2, 

Al - v; fL , A R - V% Ri (21) 

leading to the mapping: 

H MF {V df = 0) -> H MF (A = 0) 

-^eV-e d0 -e n/2 . (22) 
i<o 

This transformation, in particular, maps the FE electron 
(hole) excitations into the NB hole (electron) excitations. 

III. TRANSPORT ACROSS THE JUNCTION 

In the first subsection, we solve the BdG equations; in 
the second one, we compute the transport coefficients, in 
the limit of small electric or thermal gradient across the 
junction. 



We now find the canonical transformation that diag- 
onalizes TC MF in a self-consistent way, i.e., we solve the 
BdG equations. In the present context, the BdG ap- 
proach is essentialy the Hartree-Fock treatment of the 
two-band model. It is convenient to refer all the excita- 
tion energies to the chemical potential /i. Thus, we define 
the number operator N, 



/2 



and we replace the Hamiltonian TL MF with Ti MF ~iiN . To 
diagonalize it, we introduce the unitary transformation 



Ike 



1-kh 



= EK(j)^ + ^(j + i/2)W], 



+ w_ fc (j + l/2)/ i+ i /2 ], 



(23) 



where k is simply a quantum index, equivalent to the 
crystal momentum only in the bulk case. The idea is 
that the operators y ke , ^_ kh must diagonalize Tt MF — 
fiN and create elementary quasi-particle excitations of 
energy ui(k) (electrons) and uj(—k) (holes), respectively, 
if applied to the ground state. Therefore the equations 
of motion for the operators "fke and j—kh are 

ihj ke = [ lkei H MF -tiN]=u;(k) lke , 
ihj- kh = ["f- kh ,H MF -vN]=Lu(-k) 1 _ khl (24) 

and the Hamiltonian 7i MF — /j,N acquires the form 
H MF -^N = ]T [u>(k)4 e7ke 

k 



kh 



Note that in our definition, Eqs. (|24|) . the operator 7_ 
is a fermionic creation operator which excites a hole, 
hence the hole energy, as well as the electron energy, 
is positive, i.e. uj(k) > 0, ui{—k) > ["excitation rep- 
resentation" (ER)]. Besides, w(k), uj(—k) depend on the 
chemical potential fi and on the built-in coherence A, 
and consequently on the temperature T . The 7 opera- 
tors are analogous to the Bogoliubov-Valatin operators 
of the BCS theory; however, in the present context the 
eigenstates of the hamiltonian Tt MF have definite num- 
bers of particles, namely 7^ creates an electron and 7_jy, 
creates a hole, while the particle number of the BCS so- 
lution is an average quantity, and in that case 7^ creates 
a quasi-particle which is a mixture of electron and hole. 
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From Eqs. (|24|l and the condition that the coefficients of each independent 7-operator must be zero, we obtain 
finite-difference equations for the electron site-coefficients u and v of the junction. For the bulk semiconductor the 
BdG equations are: 

+ ft] u k (j) = -tu k (j - 1) - tu k (j + 1) - (A* - V; f ) v k (j - 1/2) - (A* + V%) v k (j + 1/2) Vj > 1, 
[u(k) +fj]v k (j + 1/2) = - (A + V <v )u k {j)- (A -V df )u k (j + 1) Vj>0. (26) 

Similarly, the hole amplitudes u and v are given by: 

[Co(-k) - p] u- k (j) = tu-k(j - 1) + tu- k (j + 1) + (A - V df ) v- k {j - 1/2) + (A + Vdf) v- k (j + 1/2) Vj > 1, 
[Q(-k) - m] v-k(j + 1/2) = (A* + V£ f ) u_ fc (j) + (A* - V#) u_ fc (j + 1) Vj > 0. (27) 



To solve Eq. I|26l) . we write down the trial two- 
component wavefunction 



Uk(j) \ = ( U k \Aka 3 

v k {j + l/2)) \v k e^ ] 



(28) 



In the following, we take without loss of generality u k > 
0, and A, V df (A L , A R , V dfL , V dfR ) real. Equation 
is compatible with l|26[) only if 



oj(k) = £ k ± E k - /i, 



(29) 



where 



a = £fc/2, 



efc = — 2t cos (fca), 



^ + |A fc -y fe | 2 , 



A fe = 2 A cos (fca/2), V fc = 2iV df sin (fca/2). (30) 

Conse- 



Since uj > 0, we choose the sign + in Eq. I|29|l 
quently, the amplitudes {u k ,v k ) are given by 

(V k -A k )u k = (£ k +E k )v k , 

plus the normalization condition 



1, 



that is 



,1 = - I 1 



I I 2 1 

\v k \ =2 



1 - 



and the phase of v k is determined by Eq. J2TJ. 
same way, the solution for holes is (ii- k > 0) 



u-k(j) 
w_ fc (i + l/2) 



W-fe 



w_ fc e- 



-ifco/2 



-ikaj 



with energy 



where 



£ fc + |A_ fc + V_*r, 



(31) 
(32) 

(33) 
In the 

(34) 

(35) 
(36) 



and amplitudes 



g-fc 
g-fc 



and the phase of determined by the equation 



(37) 



(38) 



Figure |21 shows typical quasi-particle band structures on 
both sides of the junction for the three different cases un- 
der study. In this picture the energy branches are drawn 
according to the "semiconductor representation" (SCR), 
where the quasi-particle energy e for holes is negative, 
namely 



' uj for electrons 
. — uj for holes. 



In SCR the ground state (vacuum of 7-operators in ER) 
can be represented regarding the valence band as all filled 
with electrons and the conduction band empty. Note that 
the asymmetry of conduction and valence bands close to 
the gap is remarkable in case (i) and (ii) , contrary to case 
(iii). For FE (NB) the gap is indirect, and the bottom 
(top) of conduction (valence) band is much flatter than 
the top (bottom) of valence (conduction) band, while the 
curvature of the two SC bands close to the direct gap is 
comparable. Besides, as a consequence of the transforma- 
tion 121|) . the mapping k — ► n /a — k, e — > — e transforms 
the FE electron (hole) band into the NB hole (electron) 
band. 

Once we know the semiconductor bulk energy spec- 
trum, we can also compute the probability current den- 
sity JfeeAf(j) [J-khN(J)] associated with the quasi-particle 
wavefunction (u,v) [(u, v)] (see appendix |A"|) . From 
Eq. HA11IA14II . P%|1 and we obtain 



2t 1 

JkeN(j) = —u 2 k sm(ka) + -u k \v k \ 
n n 



(A + V df ) 



x sin (fca/2 + ip) + (A - V df ) 

x sin (fca/2- tp) , v k = \v k \ e 1 ' 



(39) 
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2t 1 
J-khN(j) = —u 2 ^ k sin (ka) + -u- k \v-k\ 

x sin (ka/2 - 0) + (A - 



(A + V df ) 



x sin (ka/2 + tp) 



v_ k = \v- k \e&. (40) 



Clearly the bulk current is site-independent. It can be 
shown that 



Jkt 



N 



1 du(k) 
ah dk 



J- 



khN 



1 dto(-k) 
ah d(—k) 



( 41 ) 

i.e. the quasi-particle velocity has the same expression 
as the semi-classical one. According to formulae (|41() for 
the probability density current, the solution 1)28(1 [ 1)34) 1 ] 
represents a two-component electron (hole) wavefunction 
with wave vector k (— k) travelling from left to right if k > 
0. In the bulk we assume periodic boundary conditions 
for k and restrict its values to the first Brillouin zone, 
namely k = (2n/a) n/N s , -N s /2 < n < N s /2, N s even. 

Let us focus on case (i). In order to compute the tem- 
perature dependence of the order parameter A, recall 
definition 1 (10(1 and use the inverse transformation of 123)) 
to obtain 



A = 



^E{ fi -fcCj)e**a+i/2)(7-fch7i fch ) 

s k 

+ 40'H(i + V2)(7L7fce}}. (42) 



The average occupations (jl e ^ke/, \J-kh7-kh ) are sim- 
ply given by the Fermi function f(uj) 



(lllke) = /(*(*)) , (ll kh l- kh ) = /(«("*)) • 



with 



(43) 



= l/k B T. 



Replacing the (it, v) amplitudes with the values we have 
just computed, Eq. turns into 



A(T) 



My!^L cos{ka/2 ) 

N s 4" <■ 



[l-f(Q(-k))-f(w(k))]}, (44) 



which is the analogous of the BCS gap equation, im- 
plicitely defining A(T). The restraint d2()\i on the elec- 
tron number is an implicit definition of fJ-(T); in terms of 
7 operators it turns into: 



1 



1 



b k b k 



(45) 



We solve both Eqs. 1441) and 1)45)) simultaneously to ob- 
tain the value of A(T) and fx(T) deep inside the FE 
bulk. A critical temperature Tq exists beyond which 
A = /i = and the energy gap of the material vanishes. 



Since we know the bulk solution on the semiconductor 
side of the junction, we now study the motion of quasi- 
particles along the whole junction. The BdG equations 
for the bulk of the metal on the left side 



[u)(k) + fi] u k (j) = eVu k (j) 

-t'uktf - 1) - fUk(j + 1) , 
[u(k)+fi]vk(j + l/2)=0 Vj<0, 



(46) 



[u){-k) - fj] u-k(j) = -eVu^ k (j) 

+t'u- k (j-l)+t'u- k {j + l), 
[tS(-fc) - /*] fJ-fcO* + 1/2) = Vj<0, (47) 



have the Bloch solution 

Uk(j) 

v k (j + 1/2) 

U-k(j) 

v_ k {j + l/2) 
with energy 



(48) 



uj (q) — eV — 2t' cos (qa) — /i, 
ui'{-q) = -eV + 2t' cos (-qa) + fi, (49) 

and, from Eq. ()A9IA10)I . probability current density 

2t' 

JkeN(j) = — sin (qa), 
2t' 



J-khN(j) = — — sin(-ga), 



or, equivalently, 



Jkt 



N 



1 du'(q) 
ah dq 



J. 



khN 



1 du'(-q) 
ah d(—q) 



(50) 



(51) 



The idea is to match in some way the bulk solutions on 
both sides of the junction. To make the discussion easier, 
let us consider only the electron motion. The physical 
boundary condition is that the electron travels e.g. from 
z = — oo towards positive values of z and it is partly 
transmitted through the junction and partly reflected. 
Note that the same energy u> corresponds to two different 
bulk wave vectors q and k, while in metal/superconductor 
junctions the wave vector can be approximated by the 
Fermi vector on both sides of the interface. We always 
use k to label the coherent electronic state through the 
whole space, with the convention that both k and q, wave 
vectors in their respective bulks, correspond to the same 
energy to. It is convenient to define the wavefunctions 



*il(?) 
*ai(j + l/2) 



*2R(j + 1/2) 



' (')<• V/. (- r »2) 



= n 



Uk 



Vke 



ifea/2 



jkaj 



Vj, (53) 



8 



with the elastic scattering condition 
D (g) = Lo(k) . 



(54) 



If we dchnc 



Uk(j) 
v k (j + l/2) 



as the solution of the motion across the whole junction 
Vj, we immediately have from Eq. (|46[1 that 



Uk(j) 

v k {j+l/2) 



*2l(j + 1/2) 



Vj < 0, 



u*(0) = *il(0), 
and from Eq. if^Hjl that 



(55) 



v k (j + l/2) 



*ifl(j) 
*2fl(j + 1/2) 



Vj > 0, 



(56) 



because <|26[) and 1|46|) are linear and homogeneous. We 
have still to determine v k (l/2) and the two unknown con- 
stants T k and R k , so we need the three BdG equations 
at the interface not yet employed: 



[uo{k) + fj] u fc (0) = e dQ u k {0) - t'u k (-l) - tu k (l) - (a* l + V d ) L ) w k (l/2) . 



(57) 



+ m] »fc(l/2) = e/i/2«*(l/2) - (Al + %l) u fc (0) - (A R - u fc (l) , 



(58) 



;(fc) + M ] Ufc(l) = -tufc(0) - tu fc (2) - (Aj, - V£ fl ) w fc (l/2) - (A* + V$) « fc (3/2) . 



(59) 



We replace u fc (-l), u fc (0), u fc (l), v k (3/2), u k {2) in with *ii(-l), *il(0), * 2 b(3/2), # lfl (2), respec- 

tively, obtaining a linear system for the three unknown quantities Ufe(l/2), and T k : 



(60) 



- (A L + V dfL j R k + (u- e/i/a) «*(l/2) + (A fl - V# fl J w fc e ifcQ T fe = - (A L + 
-tfl fc + (a% - V$j?) v k (l/2) + (uu k e ika + tu k e 2[ka + (a* + V^j v k e^ ika ^j T k = -t 
(-w + e d0 - t'e^) R k + (A* + V$ £ ) Ufc (l/2) + te fc e ifcQ T fe = -u + e da - fe^'". 

To solve system (|60|) we fix uj, then we invert Eqs. (|29|l and 149(1 to obtain k and g, and hence u k and v*,: once we 
input as parameters A^, A/?, V d /L, VdfRi £do, and £/i/2, the whole coefficient matrix of system l(60j) is known and T k 
and i?fc can be eventually obtained. 



In Eq. I|52(l we chose a normalization such that the flux 
of the incident wave J]^{j) is 



jteN (i) = X ^ = \ ^ ^ < °' 

n, an, ag 

and that of the reflected wave J k c e %(j) is 

,2 2t' 



(61) 



i4r^sin(ga) < 0, (62) 
ft, 



hence, by definition, the reflection coefficient 7£(a/(g)) 
lZ(tu(k)) is given by 



I 7 refl I 
KfceiVJ 

I P nc I 
rkeN\ 



\R k \ 2 . 



(63) 



The transmission coefficient T{oS) can be calculated by 
the definition 



I Ttrans I 
\ J keN I 



'keN\ 



(64) 



with 



'keN 



! (j) = \T k \ 2 u 2 k sin (ka) + ^ |Tfc| 2 ti fc |i> fe | 



Vj > 0, (65) 



(A + V d/ )sin(fca/2 + ^) 
+ (A - V df ) sin (fca/2 - <p) 

or, more simply, by probability conservation 

T(u) = 1 - K(w) . (66) 

T(u>) is the key quantity we need to compute the ther- 
mopower. 

Because of the two-fold degeneracy of energy u), we 
have to consider also the quasi-particle associated with 
—k (and —q), with k > 0. It is clear how to rewrite 
Eq. (|52I53[) . because now the electron, coming from the 
right, is partly reflected into the right-hand side of the 
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junction and partly transmitted into the left-hand side, 
i.e. 



*2t(j + l/2) 
*2R(j + 1/2) 



U-k 



-lqaj 



(67) 



U-k 



ifea/2 



3 ifcaj 



(68) 



After that, we can proceed in the same way as above. By 
time-reversal symmetry, one has 



T(u(k)) = T(u(-k)). 



(69) 



An analogous ER procedure is used to compute the trans- 
mission and reflection coefficients T(ui) and TZ(u>) for 
holes. In this last case, Eqs. (|52I53|) turn into (with 
q,k > 0) 



*il(j) 
*2l(j + 1/2) 



R - k W e " 



(70) 



y_ kC -ika/2 



-ikaj 



(71) 



v*2flC? + V2), 
The above derivation of the transmission coefficient 
still holds for case (iii) as long as we rename and set 
the parameters as discussed in section [H] Some caution 
is needed in taking the correct sign of the wavevector k 
for certain values of parameters outside the range we ac- 
tually examined. Namely, if 2i 2 (A 2 + Vj f ) < (A 2 - Vj f ) 2 
and t 2 < A 2 — V d 2 the electron wavefunction H28|) travels 
from right to left if k > 0, like the hole wavefunction 1341) 

r2 ^ < ( A 2 



when 2£ 2 (A 2 + V 2 fl 



V 2 ) 2 ^&t 2 <V 2 -^ 2 



B. Transport coefficients and ZT 

If we apply an electric field Bora temperature gradient 
VT across the junction, we produce electric and thermal 
currents Je and Jt, respectively. In a stationary state 
the relation between fluxes and driving forces is linear, if 
fields are small enough: 37 



r Je 



(72) 



LeeE + Let^T 
Lte-E T^^VT. 

The transport coefficients L are not all independent, be- 
cause of the Onsager relation 

Lte 



L 



ET 



T 



(73) 



The coefficients are related to the electrical and thermal 
conductivities a and respectively, and to the ther- 
mopower (Seebeck coefficient) Q of the junction: 

LteLet \ 

EE ) ' 

Let 



a = L 



EE, 



Lee 



(74) 



A high value of Q alone does not necessarily imply that 
the junction is a good thermoelectric cooler, that is why 
engineers introduce the figure of merit Z: 



Z = 



Q 2 <j 



(75) 



Indeed, if we pump heat from the cooler side to the hot- 
ter one, we need high pumping efficiency (high Q), low 
production of heat through Joule heating (high er) , and 
low backwards conduction of heat (low n) & Z has units 
of inverse temperature, thus in general what is quoted is 
the dimensionless quantity ZT. We shall also use another 
natural definition of figure of merit, ZT j (I + ZT), which 
is the dimensionless figure, defined below in Eq. H93J) . 
more directly related to the transmission coefficient mo- 
menta we compute. 

Since the computation of the electrostatic and thermal 
field across the junction is complicated and not essen- 
tial to the junction thermopower, we simply assume that 
the voltage V and the temperature T are constant on 
both sides of the junction and have a sharp step at the 
interface. Thus, instead of Eq. I)72|l . we write 



Je = K EE {-SV) + K ET ST 
J T = K TE (SV) + K TT ST, 



(76) 



where we have replaced the gradients E = — W and VT 
with (SV) and ST, respectively, where SV = V (right) — 

V (left) and ST = T — T n , being T the temperature as- 
sociated to the right-hand side, and T n to the normal 
metal on the left-hand side. Since we have assumed 

V (right) = 0, we have SV = V (left) = V. The co- 
efficient Kee represents now the conductance G in place 
of the conductivity a and — (Ktt — KteKet / Kee) the 
thermal conductance in place of the thermal con- 
ductivity k, but all the other relations I|73I75|I still hold: 



K 



K 



te 



ET 



T 



(77) 



G — Kee, 

Q 



Gt — — 
Ket 



K 



KteKet\ 



K 



EE 



Q 2 G 
Gt 



(78) 



(79) 



Note that the Onsager relation (|77|l is still valid despite 
our simplifing assumptions about the V- and T-gradients. 
From Eq. I)76|l the formulae for the K transport coeffi- 
cients follow: 



E 



= ( dJ E \ 
ET \d(ST)J ST > 
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K T e = 



dJ 7 



9 (SV) J STSV=0 : 

dJ T \ 
d(ST)J ST , v=0 - 



(80) 



Here we closely follow the approach of Blonder et 
ai41 We assume that the two sides of the junction are 
in contact with perfect electron reservoirs, at different 
temperatures and chemical potentials, and that quasi- 
particle wavefunctions keep their phase coherence across 
the whole system except at z = ±oo, where they com- 
pletely loose their phase, thermalize and relax in energy 
due to inelastic scattering processes of the Fermi sea. 
Namely, we regard the contacts as perfect emitters or 
adsorbers. The electric current density, Je, will be given 
by the sum of all contributions of the quasi-particle exci- 
tations to the current, each one weighted by the correct 
Fermi distribution function, depending if quasi-particles 
originate from the left or right reservoir. Because Je is 
stationary and conserved, we can calculate it in every 
point of the space, so we choose a lattice site j < in the 
bulk on the left-hand side. In particular, the contribu- 
tion to the density current J E ~^ R produced by electrons 
going from left to right is [see Eq. 



J 



L->R 



= 2^J2 ji & C U)f("W-tV) 3<0, ( 81 ) 

s fe>0 



where the Fermi function refers to the temperature T n of 
the metal on the left-hand side, and the electrochemical 
potential differs from that on the right-hand side by eV. 
The prefactor 2 accounts for the spin degeneracy. Note 
that Eq. i|81|) gives the only electron flux running from 
left to right on the bulk metal side. On the other hand, 
the contribution to the electron current J E *~ R from right 
to left, always on the bulk metal side, is given by two 
distinct terms, one for electrons reflected at the interface 
and coming from the left side, hence in equilibrium with 
the left reservoir (with temperature T n and Fermi func- 
tion referred to eV), and one for electrons transmitted 
through the interface and coming from the right reser- 
voir at temperature T and at ground potential: 



J^ R = 2^^4f(i)/(o;(fc)- e y) 

s fc>0 

+ 2-^5>r e ans W/Mfc)) j<0. (82) 



fc<0 



Note that the second sum in Eq. I|82[l runs over negative k 
values and the related electron wavefunctions extend all 
over the junction, being therefore well defined at j < 0. 
Now we add JgTJ and using JSHl and 



k>0 



x [1 - - eV) 

+ 2^E J ^ C (i>0)T( W (fc))/(a;(fc)). (83) 



Here the notation J k n e c (j > 0) refers to electrons, incident 
on the interface, coming from z = +oo. From Eqs. 141(1 
and lf5ljl 



J 



t—*R , jL^R 
T J p. 



4e 



e duj'(q) 



N< ^— ' ah dq 

q>0 y 



x [l-n(u'{q))]f(u'{q)-eV) 

>THfc))/Mfc)), (84) 



2 1 e 9u>{k) t 



N s ^ ah dk 

fc<0 



and going to the continuum limit 
(N s a) /(27r) J dk we obtain 



J. 



L~*R | T L+-R 



2e 



7T / a 



= 2*h' dQ 



x ^MT(u'(q))f(u;>(q)-eV) 



dq 
2e 
2^h 



dk 



■k j a 



doj(k) 
dk 



T(w(k)) f(u;(k)) , (85) 



J E +J E -Xio 

x T(u>) [f(u> - eV) - f(u)] , 



(86) 



with the convention that T(u>) — if lo is not in the range 
of excitation energies allowed, and h — 2nh. Similarly, 
the hole contribution to Je is given by 



2c 



dQT(Q)[f(Q + eV)-f(Q)]. (87) 



Now we add lf%6")l and ijHIJ), passing to SCR and using the 
equivalence /(— e) = 1 — /(e), to obtain the current J E : 



2p f°° 

J E = - J de[T(e)+T(-e)] 
x [/(e - eV) - /(e)] . 



(88) 



The electron (hole) heat current density JkeT(j) 
[J-khT(J)] is defined by: 



JkeT(j) = Lu(k) J keN (j) , 
J- khT {j) = w(-fc) J-khN(j) 



(89) 



(note that in ER energies u> are referred to n). Reason- 
ing in the same way as for Je, we obtain the total heat 
current Jt'. 



2 r 00 

Jt = ^ J_de[T(e)+T(-e) 
x e [f(e-eV)-f(e)} 



(90) 



Now it is easy to compute Zf-transport coefficients from 
definitions l)80[l. We only summarize the results: 



fc<o 



2e 2 

Kee = —[—L 0l 
a 



K E t 
Ktt 



1 2e 

T T 
1 2 



Li, 



T h 



L 2 , (91) 
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with 



/oo 
de[T(e)+f(-e)] e r 
-OO 



dm 

de 



(92) 



From Eq. I|91[l one can immediately check that the On- 
sager relation l|77|) is fulfilled. Besides, note the striking 
similarity of Eq. (|91I92|) with the formalism of the semi- 
classical theory of transport (Boltzmann Equation) ^ in 
our case the role of the fc-dependent relaxation time r(k) 
is played by the transmission coefficient T, the interface 
being the scattering mechanism. From the transmission 
coefficient momenta <|92ll we compute the alternative fig- 
ure of merit 



LqLz 



ZT 
ZT+1' 



(93) 



IV. RESULTS AND DISCUSSION 

We present the results for the thermopower Q and 
the figure of merit ZT for the different types of junc- 
tion under study. Bulk parameters are chosen to de- 
scribe a large class of materials and are the same as in 
Fig. [2 with the metal band much broader than the semi- 
conductor one. In particular, for FE and NB on the 
right-hand side of the junction the indirect gap -E gap at 
T = is QAOt, i.e. one tenth of the "bare" d-bandwidth 

4* [£ gap = -t+ ^t 2 +4A 2 (0) and -t + Jt? 

respectively]. In case (iii), the gap is direct, E Eap = 
0.40i, where 4i is roughly the total bandwidth (electron 
plus hole band) of SC [E gap = 2 v /<" 2 + 4V C 2 ]. With this 
choice of parameters, the three semiconductors have the 
same gap at T = and approximately the same band- 
width. 
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A. Clean interface 

First we study the junction with a clean interface. 
With this terminology we mean that there are no d- or /- 
impurity layers (e^o = £ /i/2 = 0) at z = or z = a/2, re- 
spectively. For case (i) we assume that the change in the 
order parameter A(T) at z = is abrupt, from zero to the 
bulk value [Al — A^ ; = A(T)]. Similarly, we take the hy- 
bridization or hopping parameters at the interface equal 
to the bulk values in both cases (ii) \VdfL = VdfR — Vdf] 
and (iii) [V eL = V eR =V e ,t L = ?R = 

In Fig.^a) we plot the absolute value of thermopower 
\Q\ vs temperature. In all three cases \Q\ goes to in- 
finity as T — > 0, as expected for both indirect-gap 
narrow-band semiconductors^ and direct-gap ordinary 
semiconductors As T approaches the critical tempera- 
ture T c of the ferroelectric (k-g,T c =0.195t for our choice 
of parameters), the thermopower \Q\ of FE goes to zero, 
while \Q\ for NB and SC decreases in a exponential- like 
manner as the temperature rises. At T = Tc the FE gap 



FIG. 3: (a) Absolute value of the thermopower \Q\ vs T, for 
the three cases under study. \Q\ is in logarithmic scale. Note 
that \Q\ of FE drops to zero, as T -> T c (T c = 0.195t, with 
parameters of Fig.[5J. Inset: Magnification of the same plot, 
in linear scale, in the neighborhood of T c for FE. If T > T c 
then Q = for FE, because above the critical temperature 
the gap vanishes and the bulk FE turns into a metal [T(e) = 
T(— e)]. (b) Plot of the corresponding figure of merit ZT vs T. 
Like Q in Fig. EI a), ZT of FE goes to zero as T -> T c . Inset: 
same plot for the alternative figure of merit ZT j (ZT + 1), 
whose upper bound (ideal thermoelectric) is 1. 

vanishes and the semiconductor turns into a metal with 
symmetric bands with respect to y,, i.e. Q = 0. In the 
low temperature region, instead, for k^T < 0.05t (T < 
0.25T C ), i.e., in the region where A(T) is almost constant, 
far from the second-order ferroelectric/metal transition, 
the absolute value of Q for FE and NB is nearly iden- 
tical. For SC \Q\ is about 60%-70% of the NB value in 
the whole range of temperature. To compare our results 
with reported bulk data, we can assume a typical value 



12 




0.0 0.4 0.8 0.0 0.4 0.8 0.0 0.4 0.8 



no + n-o 

FIG. 4: Total (electron plus hole) transmission coefficient 
T(e) + T(— e) in "semiconductor representation" (SCR) as a 
function of the energy e (vertical axis). From left to right 
panels correspond to FE, NB, and SC, respectively. Curves 
for different values of the d- impurity energy Sdo are plotted. 
Note that T(e) of FE is mapped into T(-e) of NB under 
the transformation e — > — e, as a consequence of equations 
(1-H > . The d- impurity energy for four curves is indicated in 
the diagram above. Other parameters used are the same as 
in Fig. Sand T = 0. 

of At = 1 eV for FE or NB, i.e. a gap around 0.1 eV at 
T = 0, with room temperature ~ O.K. With these num- 
bers, we have \Q\ = 0.11 - 0.14 mV/K at room temper- 
ature, values comparable with those of some rare-earth 
metals with very high thermopower. If instead we as- 
sume for the broad band semiconductor SC a typical gap 
of 0.5 eV, then \Q\ is around 0.5 mV/K at room tem- 
perature, consistently with characteristic bulk data. The 
inset of Fig. |3Ja) is a magnification of the plot (in lin- 
ear scale) in the neighborhood of T c . The discontinuity 
of the derivative of Q for the FE curve is a signature of 
the second-order ferroelectric transition: if T > T c then 
Q = 0. Similar kinks are found also for the other trans- 
port coefficients G and Gt- One sees that NB and SC 
have electron transport character (Q < 0) while FE has 
hole character (Q > 0). 

To understand the above features, in Fig. 01 we plot 
the total (electron plus hole) transmission coefficient 
7~(e) +T{— e) vs energy in SCR for the three cases under 
study (at T = 0). Solid lines refer to the clean interface 
case. Since the sign and magnitude of Q is established 
by Eq. I|92|) for L\, i.e. by the competition between the 
different weights of electron and hole transmission func- 
tions T and T, respectively, it is clear that the stronger 
the electron/hole asymmetry of the transmission coeffi- 
cient, the higher the thermopower. Panels (i) and (ii) of 
Fig- El present the remarkable electron/hole asymmetry 
close to the gap, hence Q for FE has a strong hole charac- 
ter (Q > 0) while NB has a dominant electron character 
(Q < 0). The mapping of T(e) of FE into f (-e) of 



NB under the transformation e — > — e — a consequence of 
the transformation l|21|) — explains why \Q\ is the same 
for FE and NB at low T. The situation is different for 
SC [see panel (hi)], because electron and hole transmis- 
sion coefficients are more symmetric, and thus \Q\ has 
a lower value. In the limit of electron/hole symmetry 
[T(e) = T(-e), as it is the case for FE when T -> T c ], 
the thermopower is zero. 

The other transport parameters, apart from Q, are the 
electrical conductance G and the thermal conductance 
Gt- We find that, for all the three cases under study, 
there exists an activation temperature around 0.03i, be- 
low which G and Gt, as functions of T, rapidly (ex- 
ponentially) drop to zero, because the number of ther- 
mally excited carriers becomes too small. This behav- 
ior, characteristic of a material with gap, shows that at 
low temperature the main contribution to the thermal 
conductance Gt is given by the lattice, which is not in- 
cluded, thus dramatically decreasing the actual value of 
ZT. This contribution, and hence the minimum work- 
ing temperature of the junction, depends on the thermal 
impedance mismatch of the interface. 

Once all transport coefficients are known, the most 
relevant quantitity to be computed for practical appli- 
cations is the figure of merit ZT, plotted in Fig. O^b) as 
a function of the temperature. In all three cases ZT is 
a monotonic decreasing function of T. However, ZT is 
much bigger for FE and NB than for SC (but if T -> T c 
then ZT — > for FE). With the numerical parameters we 
have employed above, at room temperature ZT is ^0.3 
for FE and ~ 0.4 for NB, but already at T =100 K ZT is 
- 1.1 for FE and NB and only - 0.5 for SC. In the inset 
of Fig. EOT)), we redraw the same plot in term of the alter- 
native figure of merit ZT j [ZT + 1): this is the quantity 
we will examine in the following. While ZT has no the- 
oretical upper bound, the maximum of ZT/{ZT + 1) is 
1, corresponding to ZT = oo. 

It is likely that the one-dimensional model artificially 
enhances ZT . Besides, we do not know how the order 
parameter A(T) actually varies at the interface, and all 
the effects of charge polarization are neglected: this does 
not necessarily imply a reduction of ZT . The only scat- 
tering mechanism we have considered is the interface: in 
particular, the contribution of phonons to the thermal 
conductance is neglected. 



B. d-impurity layer 

Now we study the effect of an overlayer of d-impurity 
atoms at the interface. This structure could be built by 
epitaxial growth techniques: we find that this configu- 
ration does not improve considerably the figure of merit 
ZT. 

In our one dimensional virtual crystal model, we de- 
scribe for simplicity the overlayer by putting one atom 
at z — with site energy Sdo- Since in the calculation 
we leave unchanged the hopping coefficient t and all the 



13 



other parameters with respect to the clean-interface case, 
we also require, from a physical point of view, that the 
atomic orbital of the impurity is still of d-type. The es- 
sential point here is that the atom at z — substantially 
participates in the electronic motion, contrary to a /-site. 
In case (iii) this distinction between /- and rf-sites is no 
longer relevant: actually we look at SC only for compar- 
ison. One could also think of growing an overlayer at 
z = only partially filled with impurities: in that case 
the values of transport coefficients should be an inter- 
polation between the two limiting values corresponding 
to the cases of 0% (clean surface) and 100% impurity 
concentration within the layer. 

To understand the role of the d-impurity layer in the 
transport, consider the transmission in Fig. Here, as 
well as the total transmission coefficient vs energy for the 
clean interface (solid lines), we have also plotted curves 
corresponding to increasing values of the impurity level 
Edo- As the impurity level energy e^o increases, the trans- 
mission is uniformly depressed over the whole range of 
energies. In case (i) and (ii), results depend only on the 
absolute value of £do, while (iii) is more complex, due to 
the presence of a second nearest neighbor hopping coef- 
ficient. 

While the depression effect of the impurity seems im- 
portant for both NB and FE semiconductors, some cau- 
tion is needed in accepting these results. First, one can- 
not make le^ol arbitrarily large, because its value is physi- 
cally limited and cannot differ too much from typical bulk 
energies, otherwise the impurity would be screened. Sec- 
ond, in our computation we have taken all parameters but 
£do unchanged with respect to the clean-interface case, 
and it is clear that this approximation becomes worse as 
long as \edo\ increases. For example, A/ or Vdjh would 
surely change as e^o varies. 

The conclusion is that the enhancement of the figure of 
merit ZT j (ZT + 1) for reasonable values of e^o is quite 
limited. In general, as long as we increase the magnitude 
of Edo, we slightly enhance ZT uniformly over the whole 
range of temperatures. It is remarkable that FE and NB 
values of ZT are always much higher than SC values. For 
example, at room temperature, with the usual choice of 
numerical parameters and Sdo — 5i, ZT ~ 0.4 for FE, 
-0.6 for NB, and -0.2 for SC, but already at 100 K 
ZT ~ 1.8 for FE and NB, while it is only - 0.4 for SC. 



C. /-impurity layer 

In order to dramatically improve the figure of merit, we 
propose the insertion of a rare-earth overlayer at z = a/2, 
with /-level energy e/1/2. Contrary to the d- impurity 
layer discussed in the previous section, here it is essential 
that the localized impurity level is of /-type, i.e. hardly 
sharing the electronic conduction. This is not the case 
for SC, where tunneling from this "/-site" to adjacent 
neighbors is allowed: in fact, for SC, the situation is ba- 
sically similar to the previous (i-impurity case. In case 
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FIG. 5: Total (electron plus hole) transmission coefficient 
T(e) + T(— e) in SCR as a function of the energy e (vertical 
axis). From left to right vertical panels correspond to FE, 
NB, and SC, respectively. Curves for different values of the /- 
impurity energy e f 1 / 2 are plotted. A consequence of equations 
12 It is that T(e) of FE for e/1/2 > can be mapped into 
T( — e) of NB under the transformation — e — * e, — e/1/2 — * 
e fi/2- The parameters used are the same as in Fig. with 
T = 0. 



(iii), hopping to adjacent neighbors is allowed both from 
the "/-site" and from the "d"-site. In fact, SC results for 
the figure of merit ZT are not dissimilar from the values 
obtained in the d-impurity case. In order to gain some 
insight into the thermopower behavior, in Fig. [3] we have 
plotted the total transmission coefficient 7~(e) + T(— e) 
vs e for different values of £/i/2) similar to Fig. 0] FE 
and NB curves [panel (i) and (ii), respectively] are quali- 
tatively different from those of SC [panel (iii)], as we set 
£fi/2 to negative values in the energy band range. For FE 
and NB, T goes to zero in the neighborhood of e/1/2, as 
if the hole were completely backscattcrcd from the inter- 
face when resonates the energy with that of the impurity 
atom. On the contrary, in case (iii) the effect is opposite, 
with T gaining weight for energies close to e/1/2 ■ The 
trend is similar for £/i/2 > 0, showing total reflection 
in the neighborhood of positive values of e/1/2 1 m case 
(i) and (ii). These results demonstrate that while the 
/-impurity level of FE and NB does not share the elec- 
tronic conduction, due to localization, it strongly affects 
the transmission, because either in case (i) the quasi- 
particle excitation is a coherent superposition of d- and 
/-states, or in case (ii) the "bare" bands are hybridized. 
The case of SC is excluded here. 

The overall effect of the /-impurity layer on transport 
is so strong that it even changes the dominant (electron 
or hole) character of the thermopower, i.e. the sign of Q. 
For the sake of simplicity, we now consider only FE. In 
Fig. we plot Q vs T for different values of £/i/2- Left 
panel refers to e/1/2 < 0. We see that, as we set e/1/2 
from zero to negative values, Q changes sign: already at 
e/1/2 = -0.lt Q has electron character (Q < 0) as T — > 0, 
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FIG. 6: Thermopower Q vs T for FE. Left panel: Curves 
for different values of the /-impurity energy £/i/2 (negative). 
Right panel: The same, with e/1/2 > 0. Note that Q changes 
sign for some curves as T varies, and that \Q\ — > 00 as T — > 0. 
Bulk parameters as in Fig. [5] 



while for fc B T > 0.02i Q has hole character (Q > 0). 
This behavior can be understood by examining Fig. [3] 
the impurity level drastically diminishes the weight of the 
transmission coefficient at the top of the valence band, 
while increasing it at the bottom of the conduction band, 
so that the sign of Q is reversed at low temperatures, 
where the only excited carriers are those whose energies 
are close to the gap. Note that Q is extremely sensitive 
to the position of e/1/2, m contrast to the case of the 
d-impurity layer. The right panel of Fig. presents a 
similar situation for e/1/2 > 0. Here the situation is a 
little less obvious: as we raise the value of e/1/2 first Q 
rises then drops and changes sign. This is due to the 
asymmetry of electron and hole bands. At first e/1/2, 
now positive, is close to the bottom of the conduction 
band, favoring the hole transport because it depresses 
the thermally activated electronic channels. Then, as 
e/1/2 is increased, weight is added to the transmission 
coefficient at the bottom of the conduction band. This 
weight is "swept away" from the energy neighborhood 
resonant with e/1/2, that now is higher in energy with 
respect to the band bottom. This mechanism also makes 
Q change sign. 

To study the behavior of ZT as the /-impurity 
e/1/2 varies, we take a series of "snapshots" depicting 
ZT I {ZT + 1) vs T for different values of e /1/2 in Fig. □ 
Figure EJa) focuses on e/1/2 < 0, in a range between 
and -0.2t. We see that, as e/1/2 is lowered from its 
zero value, ZT/(ZT+l) conspicuously decreases: in 
this range Q changes sign. At e/1/2 = -0-1 and -0.125i 
a temperature exists at which ZT = 0: these tem- 
peratures correspond to the zero of Q in Fig. HO (left 



panel). At e /1/2 = -0.2t, however, ZT / {ZT + I) be- 
gins to rise. In Fig. Efb) we analyze the "high-value 
region" of ZT/{ZT+1). At e/i/ 2 = -0.24i the curve 
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FIG. 7: (a) ZT/ {ZT + 1) vs temperature T for different 
values of the /-impurity level energy e/1/2 for FE. Bulk pa- 
rameters as in Fig. (b) The same as Fig. EI a ), f° r lower 



values of e 



/1/2- 



Note that the "shoulder" in the curves is 



shifted towards higher values of T, as e/1/2 is decreased. The 
dashed line represents the envelope of the curves. 



is quite depressed around /cbT =0.1 (room temperature 
with the usual numerical parameters), but it acquires gi- 
ant values around T = 0, close to the theoretical upper 
bound 1 . Moreover, the shape of the curve is very differ- 
ent from the typical pattern of the c?-impurity layer: in 
this latter case the curve is convex around T = 0, while 
in the present situation ZT j {ZT + 1) is concave, i.e. if 
one slightly departs from T = ZT will still be very 
high. This approximate flatness of the curve suggests 
that a stable low-temperature working point should ex- 
ist for a junction-based device. As e/1/2 is furtherly de- 
creased, ZT I {ZT + 1) shows a concave "shoulder" which 
is shifted towards higher temperatures. This feature im- 
plies that, once we fix a temperature, it is possible to find 
an optimal value of e/1/2 maximizing the figure of merit. 
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FIG. 8: Maximum value of ZT / (ZT + 1) vs T for an opti- 
mized choice of the /-impurity level energy £/i/2- Case (i) 
(solid lines) and (ii) (dot-dashed lines) are represented. The 
up (down) triangles refer to the positive (negative) values of 
local maxima £/i/2, computed at fixed T. Note that the solid 
line with down triangles is the envelope curve of Fig. EJb). 
Bulk parameters as in Fig. [5] 



This locus of optimal working points is manifestly given 
by the envelope curve for the "shoulders" (dashed line). 
The behavior of ZT j (ZT + 1) as £/i/2 rises from zero 
( e /i/2 > 0) is similar, but the sequence is inverted with 
respect to the case we have just examined: this time first 
ZT j (ZT + 1) dramatically increases, and then rapidly 
drops to low values. 

Figure [H] shows the maximum values of ZT j (ZT + 1) 
vs T that can be reached with a suitable /-impurity en- 
ergy £/i/2- Up (down) triangles indicate positive (neg- 
ative) optimal energies £/i/2- Analogous results for NB 
(dot-dashed lines) are shown. The curves correspond to 
the loci of the T-£/i/2 space like the envelope curve of 
Fig. [7fb). From Fig.|S]we can derive the absolute max- 
ima of ZT. Our results are that, at room temperature 
(ksT = O.li, with the usual numerical parameters) ZT ~ 
1 for FE and 1.6 for NB, at T = 150 K ZT is - 6 in both 
cases, at T = 100 K is - 13, and at T = 40 K is - 100. 

These extremly high values suggest the possibility of 
engineering a periodic lattice of (5-layers made of metal 
and strongly correlated semiconductor to create an effi- 
cient thermoelectric device. Borrowing an intuition of 
Refs. |lj and [lja, the bias applied to the superlattice 
should be such that the electronic transport occurs per- 
pendicular to interfaces. The thickness of each layer 
should be comparable to the electronic mean free path, 
and large enough to prohibit electrons from tunnelling. 
In this regime the electronic motion within each layer 
is ballistic, and the interfaces can be considered as the 
only scattering mechanisms. Mahan and Woodsi^ call 
thermionic a similar device with an ordinary semicon- 



ductor replacing the strongly correlated one, because the 
electronic motion is ballistic and the expression for the 
current is the Richardson's equation. However, here the 
situation is different, because the transmission coefficient 
has such a sharp variation with energy that the current 
has a more complex form than the Richardson's expres- 
sion and it requires a full analysis of the role of the in- 
terface in transport. The idea here is to grow a /-atom 
overlayer at z = a/2 with optimal /-level energy £/i/2 
depending on the working temperature of the device. In 
addition, we find that the magnitude of thermal conduc- 
tance changes only slightly with £/i/2- The <5-layer of 
rare-earth atoms could also act like an additional source 
of strong scattering for phonons, thus lowering the lat- 
tice thermal conductance.- Therefore the junction is a 
promising candidate as a thermoelectric device. 

D. Effects of the surface on the lattice structure 

We briefly comment on the influence which the inter- 
face has on the ideal lattice structure. We consider the 
narrow band semiconductor NB. If the effect of the sur- 
face at z = is to locally shorten the lattice constant 
a, this will certainly change the hybridization parame- 
ter VdfL between the d-site at z = and the /-site at 
z = a/2. Presumably VdfL will be increased. Modifying 
nothing but Vdf l with respect to the clean interface case, 
we find that ZT is only slightly affected by this relaxation 
effect. In particular, to change VdfL U P to 10% almost 
rigidly shifts ZT/ (ZT + 1) over the whole temperature 
range by approximately 10 - 15 %. From a check of the 
transmission coefficient, we find that this surface effect 
represents only a minor perturbation to the electronic 
transport across the junction. 

V. CONCLUSIONS 

We have made a qualitative theoretical study of the 
possibility that a junction of metal and FE or NB (as op- 
posed to bulk materials or to a junction metal/SC) pro- 
duces high thermopower. This is possible if a Mayer of 
suitable rare-earth impurity atoms substitutes the origi- 
nal layer at the interface. The localized character of the 
impurity /-orbital has a strong effect on the transmission 
of carriers across the junction. The figure of merit ZT 
attained is very high, especially at low temperatures. A 
realistic device is proposed which exploits the thermo- 
electric potentialities of the junction. 
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APPENDIX A: CURRENT OPERATOR 

In this appendix we derive the form of the density 
current operator J for the effective Hamiltonian Tt MF 
of Eq. I|18fl (or 7i MF — /i). This is a key quantity in com- 
puting the transport properties of the junction. 

One way to find J is to exploit the gauge invariance of 
the hamiltonian 118|) . In the presence of an oscillating, 
uniform electric field E(t) = Ee~ luJt along z (t is the time 
and ui/2tt is the frequency), one must add to H MF a term 
representing the coupling to the applied field. One thus 
introduces the vector potential A(t) = cEe~ luJt /iui such 
that E(t) = —A/c, where c is the speed of light. The 
gauge transformation 



dj = dj exp 
/j + l/2 = /j + i/2exp 



-A(t)aj 



he 



A(t)a(j + 1/2) 



, (Al) 



removes the coupling to the vector potential A(t) pre- 
serving the gauge invariance of Tt MF ; if we apply the 
transformation (|A1|) to 7Y MF we obtain a new Hamilto- 
nian H. MF ' which must be given, expanding to first order 
in A, by 



n MF ' ~ H MF - —A. 

c 



(A2) 



Equation i|A2|) permits to identify J, whose expression is 
•/ = E J 0')' ( A 3) 



j(j) = -et'd^d, +H.c. Vj<0, (A4) 



Here J(j) is the electric current operator at atomic site j 
(in one dimension the current coincides with its density). 
If the transport across the junction is stationary, at each 
site j the current must be conserved; thus, to identify 
the density current associated with each quasi-particle 
excitation, we can calculate the operator J(j) at some 
suitable atomic site j where its expression is simpler: this 
is used in Sec. IIIII For j < 7i MF is a simple tight- 
binding Hamiltonian, and the formula for J(j) becomes, 
replacing d and / operators in l|A4|l with 7 operators: 



2et' 



^) = ^E Im [^>^' + 1 )wL 



Ike 



2et' 



— E lm [u-k(j) u*_ k (j + 1)] ji kh i-kh 



2et' 



+ ^E Im ^0')^(j + 1 )]- ( A7 ) 



It turns out that the constant third term on the right- 
hand side of Eq. (|A7|I is always zero, hence it is easy to 
identify the electric current Jke(j) and J-kh{j) associ- 
ated with each electron or hole excitation, respectively: 



j co = ^E[ J ^'H< 



Ike 



J-kh(j)ll kh 1-kh , 



(A8) 



2et' 

i(j) = T ImK(j)%(j + l)] Vj<0, (A9) 



2et' r_ 
J-kh(j) = -— Im u-k{j) 

x u* fc (i + l)] Vj<0. (A10) 

The particle current density for electrons JkeN(j) is given 
by 



JkeN (j) 

and for holes J-khN(j) by 



n 1 



J. 



khN 



(J) 



Jke(j) 



J-kh(j) 



(All) 



(A12) 



j ( o) = r td[d + 2n 



(A L + V d f L ) f! /2 d 



Afl-V^fl) fl /3 dx +H.c, (A5) 



J(j) = ~et4+i d J + ^ [ ( A + v «) f]+i/2 d i 

- (A - V df ) /] +1/2 di+il + H.c. Vj > 0. (A6) 



The picture consistent with these results is that 
(uk(j) , Vk(j + 1/2)) [or ("S,w)] represents the two- 
component electron (hole) wavefunction of the elemen- 
tary excitation, and JkeN(j) (J-khN) is its probability 
current density. Indeed, equations l|A9IA10l) can be de- 
rived also from the continuity equation of wavefunctions, 
if one associates the electron amplitude (u, v) with the 
charge e, and the hole amplitude (u,v) with — e. To 
obtain the continuity equation, one has to multiply the 
time-dependent BdG equation gSJ) [l(17|l] by u* k (j,t) [«*], 



and the complex conjugate of (|46|l [H47JI] by Uk{j,t) [u], 
and finally subtract one term from the other. 

In an analogous way we can also write the expression of 
the current for j > 0. If A and Vdf are real, one obtains: 

let e 
Jke(j) = -j-hn [ul(j) u k (j + 1)] + - (A + V df ) 

x ImK(j> fe (i + l/2)]-|(A-%) 

x ImK(j + l)^(j + l/2)], (A13) 

let e 
J-kh(j) = — jplm [u-k(j) u*_ k (j + 1)] - - (A + 

x Im [u. k (j) v*_ k (j + 1/2)] +|(A- 

x Im[u_ fe (i + l)C fe (i + l/2)]. (A14) 
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